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Based on recent works [1,2], we present the results of calculations for sev- 
eral physical quantities (meson masses, the S parameter, etc.) in a vectorial 
gauge theory, as a function of the number of fermions, Nf. Solutions of the 
Schwinger-Dyson and the Bethe-Salpeter equations with the improved ladder 
approximation are used for the calculations. We focus on how the values of 
physical quantities change as one moves from the QCD-like (non-walking) to 
walking regimes. 



1. Introduction 

We consider a (3 + l)-dimensional vectorial gauge theory (at zero tempera- 
ture and chemical potential) with the gauge group SU(iV c ) and Nf massless 
fermions transforming according to the fundamental representation of this 
group. For N c — 3, if one took Nf = 2, this would be an approximation to 
actual QCD. We restrict here to the range Nf < (ll/2)N c for which the 
theory is asymptotically free. An analysis using the two-loop beta function 
and Schwinger-Dyson equation leads to the inference that for Nf in this 
range, the theory includes two phases: (i) for < Nf < Nf :Cr a phase with 
confinement and spontaneous chiral symmetry breaking (SxSB); and (ii) 
for Nf^ cr < Nf < (11/2) N c a phase with no spontaneous chiral symmetry 
breaking (plausibly a non-Abelian Coulomb phase). We shall refer to Nf^ cr , 
the critical value of Nf, as the boundary between these two phases. 

For Nf slightly less than Nf^ cr , the theory exhibits an approximate 
infrared (IR) fixed point with resultant walking behavior. That is, as the 
energy scale \x decreases from large values, a = g 2 / (4n) (g being the SU(AT C ) 
gauge coupling) grows to be 0(1) at a scale A, but increases only rather 
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slowly as fi decreases below this scale, so that there is an extended interval 
in energy below A where a is large, but slowly varying. Associated with 
this slowly running behavior, the resultant dynamically generated fermion 
mass, E, is much smaller than A. In addition to its intrinsic field-theoretic 
interest, this walking behavior has played an important role in theories of 
dynamical electroweak symmetry breaking [3]- [9]. As Nf approaches Nf^ cr 
from below, quantities with dimensions of mass vanish continuously; i.e., the 
chiral phase transition separating phases (i) and (ii) is continuous. Recently, 
meson masses and other quantities such as the generalized pseudoscalar 
decay constant fp and the S parameter [10] were calculated in the walking 
limit of an SU(AT C ) gauge theory [11,12]. 

It is of interest to investigate how meson masses and other quantities 
change as one decreases Nf below Nf tCr , moving away from the boundary, 
as a function of Nf, between phases (i) and (ii), deeper into the confined 
phase. For this purpose, in Ref. [1], as in Rcfs. [11,12], we use the Schwinger- 
Dyson (SD) equation to compute the dynamical fermion mass E (general- 
ized constituent quark mass) and then insert this into the Bethe-Salpeter 
(BS) equation to obtain the masses of the low- lying mesons and other quan- 
tities. We restrict to an interval of Nf values for which the theory has an 
infrared fixed point. For dcfinitencss, we take N c — 3; however, iV c enters 
only indirectly, via the dependence of the value of the infrared fixed point 
a* on N c . Hence, our findings may also be applied in a straightforward way, 
with appropriate changes in the value of a*, to an SU(AT C ) gauge theory 
with a different value of N c . 

In order to study meson masses and other quantities as one moves away 
from the boundary between phases (i) and (ii) , it is first necessary to know 
as accurately as possible where this boundary lies, as a function of Nf, 
i.e., to know the value of Nf^ cr . For sufficiently large Nf, the beta function 
(calculated to the maximal scheme- independent order, namely two loops) 
has an IR fixed point at 

-4tt(117V c - 2N f ) 
34 AT? - 13N c N f + 3Nc L N f 

Requiring that a* be sufficiently large as to yield spontaneous symmetry 
breaking in the context of an approximate solution to the SD equation for 
a fermion yields the condition that Nf < Nf^ cr , where [9] 

2iV c (5(W e 2 - 33) 
' V/ ' cr " 5(57V c 2 -3) • [) 

For N c = 3 this gives Nf iCr ~ 11.9. These estimates are only rough, in view 
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Fig. 1. Numerical solutions for S (left panel) and fp (right panel), for several values of 
a* (indicated by ()). For comparison, we show eq. (3) with c = 4.0 (left panel), c = 1.5 
(right panel) from a fit to the results for 0.89 < a* < 1.0. 



of the strongly coupled nature of the physics. Effects of higher-order gluon 
exchanges and instantons have been studied in Refs. [13]. 

In our analysis, what we actually vary is the value of the approximate IR 
fixed point a*, which depends parametrically on Nf. Thus, although our SD 
and BS equations are semi-perturbative, the analysis is self-consistent in the 
sense that our a cr really is the value at which, in our approximation, one 
passes from the confinement phase with S%SB to the chirally symmetric 
phase, and our values of a» do span the interval over which there is a 
crossover from walking to QCD-likc (i.e., non- walking) behavior. 



2. Schwinger-Dyson Equation 

We first use the Schwinger-Dyson equation for the fermion propagator to 
calculate the dynamically generated mass £ of this fermion. In Fig. 1 (left 
panel) we show the solution for the dynamical fermion mass £ as a function 
of a». A fit to the numerical solution in the walking region 0.89 < a* < 1.0 
[11] found agreement with the functional form 

E = cAcxp -7r(— -l) , (3) 

L ^ttcr ' 

with c = 4.0 (see also Refs. [6,9]). Our calculations for larger ct* show the 
expected shift away from walking behavior. This shift is evident in Fig. 1 
for a* larger than about 1.2. Our calculation of S, shown in Fig. 1, shows 
that S/A increases substantially, by about a factor of 30, from a value of 
about 0.01 at a* = 1.0 to 0.32 at a* = 2.5, much closer to the value of 0(1) 
for this ratio in QCD. 

Another quantity of interest is the pseudoscalar decay constant fp, the 
Af/-flavor generalization of the pion decay constant, f n . In Fig. 1 (right 
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Fig. 2. Values of meson masses divided by A calculated from the Schwingcr-Dyson and 
Bethe-Salpeter equations in the range 0.9 < ct» < 2.5. 

panel) we show our results for fp calculated by substituting our solution 
for E(/c 2 ) into the Pagels-Stokar formula. In the walking limit, fp has been 
shown to satisfy a relation similar to eq. (3), i.e., it is exponentially smaller 
than the scale A. We display, as the dotted curve, the fit from Ref. [11] for 
the walking interval 0.89 < a* < 1.0, given by eq. (3) with c = 1.5. Our 
results show the change from this walking type of behavior as a* increases 
above this range; specifically, as a* increases from 1.0 to 2.5, fp/A increases 
substantially, from about 3 x 10~ 3 to about 0.08. This is similar to the factor 
by which we found that E/A increased as a* increased through this interval. 

3. Calculation of Meson Masses 

We next present the results of the numerical calculations for meson masses, 
obtained by solving the homogeneous BS equation [1]. (As in Ref. [11], 
we have checked and confirmed that the flavor-adjoint pseudoscalar meson 
mass is zero to within the numerical accuracy of our calculation.) In Fig. 2, 
we show the values of meson masses divided by A calculated from the SD 
and BS equations. In Fig. 3 we plot the values of Ma,v,sI fp (left panel) and 
Ma,s/Mv (right panel). Here, the subscripts V, A and S represent vector, 
axial-vector and scalar, respectively. 

Our calculations yield a number of interesting results. We summarize 
these for the changes in these meson masses as a* increases from 0.9 to 
2.5. The ratios of the meson masses divided by A increase dramatically, by 
factors of order 10 2 , approaching values of order unity at a* = 2.5. This 
amounts to the removal of the exponential suppression of these masses 
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Fig. 3. Values of meson masses divided by fp (left panel) and Ma,s/Mv (right panel) 
calculated from the Schwingcr-Dyson and Bcthc-Salpcter equations. 

which had described the walking limit near Nf tCr , as one moves away from 
this limit into the interior of the confined phase. For example, Mg/fp 
increases monotonically from about 4 to 7, thereby approaching to within 
about 35 % of the value 10.7 in QCD for M ao /f n , while My/fp decreases 
from about 11 to 9, rather close to the value 8.5 for M p j / w and M u j 'f n in 
QCD. The ratios M A /M V and M s /M v , which were found in Ref. [11] to 
have values close to 1.0 and 0.36, respectively, in the walking limit, both 
increase in the interval of a* that we study reaching about 1.2 and 0.74, 
respectively, at a* = 2.5. For comparison, these ratios are approximately 
1.6 and 1.3 in QCD. 



4. Calculation of the S parameter 

In this section we present the results of our calculations of S [2] in the 
crossover region of the theory between the walking limit at a* \ a cr and 
larger values of a* that move toward the QCD-like regime. Here, S rep- 
resents the contribution to the S parameter from one fermion isodoublet. 
(Studies of S in the walking limit include [12,14].) We calculate S via the 
relation S = 4ir(IL vv (0) - W AA (0)), where Il vv (q 2 ) and U AA (q 2 ) arc the 
vector and the axial-vector current-current correlation functions. These cor- 
relators are computed by solving the SD equation and the inhomogeneous 
BS equation [2]. 

In Fig. 4, as a function of a*, we plot the value of S n , the value of S 
normalized by its value at a* = 1.8, namely, 0.47. This figure shows that 
S n , and hence also S, decreases by about 40 % as a* is reduced from 1.8 
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Fig. 4. Plot of S n for several values of a* in the range of 0.9 < a* < 1.8. As indicated 
by the subscript rt, the values are normalized by the value of S at «* = 1.8, i.e., 0.47. 

to 0.9, or equivalently as Nf is increased from 10.3 to 11.6. Reinserting the 
factor of the number of fermion isodoublets, Np = Nf/2, to get S itself, we 
obtain a decrease by about 30 % in S, since Nd only increases by about 10 
% over this range. Thus, our calculation shows that for this range of values, 
S decreases significantly as one moves from the QCD-like to the walking 
regimes. We recall that the (improved) ladder approximation to the BS 
equation can overestimate S in QCD by as much as 30 % [15]. Hence, in 
addition to the demonstrated decreasing trend of S and S as a* decreases 
from 1.8 to 0.9, one may, separately, comment that the absolute magnitude 
of these quantities could be about 30 % smaller than the values yielded 
by our ladder approximation. Our results thus strengthen the evidence for 
the reduction of S in a walking, as opposed to QCD-like, gauge theory, 
and are relevant to assessing the impact of the S'-parameter constraint on 
technicolor theories. 

5. Summary 

In summary, using numerical solutions of the Schwinger-Dyson and Bethc- 
Salpeter equations, we have calculated several physical quantities, including 
fp, meson masses, and the S parameter, as a function of the approximate 
infrared fixed point, a*, or equivalently, the number of massless fermions, 
Nf, in a vectorial, confining SU(7V) gauge theory. Our results show the 
crossover between walking and non-walking behavior in a gauge theory, 
and demonstrate that S and also S decrease significantly as a* decreases 
in this range. 
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